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Within the parity doublet model coupled to the linear sigma model including vector mesons it is 
possible to describe vacuum properties of the low energy mesons and to achieve nuclear saturation 
at nonzero density. Motivated by recent studies we investigate the emergence of inhomogeneous 
condensation in the parity doublet model. As a first step the chiral density wave (CDW) is con- 
sidered, which allows for a straightforward investigation of inhomogeneous condensation. As a 
result it is still possible to have a homogeneous ground state of nuclear matter, but at larger baryon 
chemical potential the CDW is favored with respect to the homogeneous phase. 
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1. Introduction 

For many decades the phase diagram of the quantum chromodynamics (QCD) was believed to 
be simple |Jlp. Only two phases were considered: a confined phase, where the relevant degrees of 
freedom are hadrons which are build from quarks and gluons, and a deconfined phase, where the 
relevant degrees of freedom are quarks and gluons [Q. 

With increasing computational power, lattice simulations were able to test the above mentioned 
picture at nonzero temperature and zero quark chemical potential [|3|, Q]. Since, at present lattice 
simulations are not able to access the high quark potential regime, effective models and QCD-like 
theories have to be used. These approaches suggest a different picture, according to which a rich 
phase diagram with a complicated structure is realized. At low temperatures and at low densities 
matter is confined and the relevant degrees of freedom are hadrons; on the contrary, at very high 
densities a color superconducting phase should be present [Q], but it is unclear how the transition 
from nuclear matter to such phases looks like. Also the critical quark chemical potential for a 
transition to such a phase is unknown. Beside these uncertainties, there is also the possibility that 
further states of matter e.g. quarkyonic matter [^J exist. 

All approaches to describe QCD at finite density have to use simplifications. For example 
the NJL model [Q] takes only quark degrees of freedom into account and therefore is not the best 
choice to describe the regime of confined matter. On the other hand, models based on hadrons have 
no chances to access the regime where quarks and gluons are the relevant degrees of freedom. 

The general idea of inhomogeneous condensation goes back to Migdal [^]. He first introduced 
inhomogeneous condensation within nuclear matter which was realized via the famous chiral den- 
sity wave (CDW). The idea of inhomogeneous condensation never got out of fashion and only a few 
years ago the Gross-Neveu model in two dimensions could be solved exactly, analytically []|] and 
numerically It has be shown that for high densities and low temperatures inhomogeneous con- 
densation dominates the phase diagram. Based on these results it was shown that the NJL-model 
in four dimensions shows the same class of modulations as the Gross-Neveu model [jKj. 

It is of major interest to also look for the possibility of inhomogeneous condensation within a 
fully hadronized model that is able to describe vacuum properties and to describe nuclear matter 



saturation. The parity doublet model UllQ successfully reproduces vacuum phenomenology [ |12| ] 
and physics at finite density [13]. The next step, achieved here, is to test the ground state at nonzero 
density for the formation of a inhomogeneous condensation. Furthermore a generalization of the 
parity doublet model is straightforward and also exotic matter states like tetraquarks and glueballs 



can be introduced []14[]. In this work however, we use the model studied in Ref. [J13p . 
2. The Model 

The mesonic part of the Lagrangian is a SU (2) linear sigma model including vector mesons 



[15]. It has the following form: 



1 1 A 

+ -m 2 {o 2 + n 2 ) + -m 2 (0 (O i;L (0^ - -(o 2 + K 2 ) 2 + EO , (2.1) 
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with the field strength tensor for the vector meson fields F^ v = d^CDy — d v co^. In the vacuum chiral 
symmetry is spontaneous broken; this is achieved in the model via the sign of the parameter nr. 
The vacuum exception value (v.e.v.) of the field a is nonzero and its value Co = <p corresponds at 
zero temperature and density to the pion decay constant, q> = f n . 



The baryon part of model considers, besides the nucleon N, also its chiral partner N* [|13p. The 
Lagrangian is formulated in terms of the bare fields and ife: These fields are chiral eigenstates 
but not mass eigenstates of the model. The physical fields N and N* emerge as superpositions of 
\f/i and \jf2- The fields Xj/\ and 1//2 transform under SUi{2) x SU (2)r the following way: 

¥i,r ->• U R ¥\,r , U L ¥\,l , ¥2,r -> U L Y2,R , ¥2,l -> U R y 2 ,L • (2.2) 

Notice that \\f 2 transforms in a mirror way with respect to the field Yi- Besides the well known 
kinetic terms and the coupling to mesons, the mirror assignment allows to construct a further chiral 
invariant bilinear term: 

Y2,L ¥\ ,R ~ ¥2,R ¥l X ~ ¥l ,L ¥2,R + ¥l ,R ¥2,L = ¥2 75 ¥l ~ % 75 ¥2 ■ (2-3) 

The full Lagrangian in the parity doublet model is: 

Jzf B = yWy/i - -g\ ¥1 {o + 275? • n) ¥\ + yWy>2 - -jgify ( a - *7s? ■ ^) ¥2 

-gw } ¥imoh¥i -g { a¥2i7oOX)¥2 + mo(¥275¥i ~ ¥i75¥i) + -&m ■ (2.4) 

The term proportional to generates an additional contribution to the mass of the nucleons as 
well as a mixing between Yi an d ¥2- This leads to the aforementioned difference between the 
chiral and mass eigenstates. In fact in order to obtain the physical masses of the nucleons, ¥\ an d 
\ff 2 have to be transformed according to the following transformations: 




1 / exp(5/2) 7 5 exp(-5/2) 




V2cosh<5 \ v 7sexp(-5/2) -exp(<5/2) 
Resulting from this transformation the now physical masses of the Af* and A^ resonances read: 

1 In ~ (p 2 +^i+U^g i -^g 2 ]( P 



(2.5) 




™ N * = ^\l[^gi + ^§2 ) <P 2 +*m( ) - l - Qgj - ^g 2 J <p . (2.6) 

From the equations for the nucleon masses, the relevance of the parameter mo is clear. It allows 
for a mass of the nucleon even if the chiral symmetiy is restored, (p = 0. The chiral condensate 
is not solely responsible for the mass of the baryons, but generates the mass splitting between the 
nucleon and its chiral partner. Sending the parameter — > one reobtains the naive assignment, 
and the masses are generated only by the chiral condensate. 

The parity doublet model combines the linear sigma model with the nucleon Af and it's chiral 



partner Af*. In the mean field approximation [16] the Lagrangian reduces to the potential: 



1 1 A 

f M = --mV--m^ 2 + -9 4 -8ip 1 (2.7) 
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and the grand canonical potential in the no see approximation reads: 

^/V = y M + J £j^y fj 3 k{Ej(k)-^)d{Ej(k)-^) , (2.8) 

where the sum runs over the resonances N and N* with the corresponding fermionic degeneracy 



factor Yj. The energy has the known form Ej = Jk 2 + mj. The 6 function requires the chemical 



potential to be jtij = ju ; - — g w (Oo = yk 2 + m2 j- Minimizing £1 with respect to (p and Oq leads to the 
nonzero density behavior. This model has been studied and extended in Refs. [|l2|, 14]. 

Following the promising results in the above mentioned publications, the model has to be 
tested for inhomogeneous condensation. A straightforward approach is the CDW, which can be 
parameterized in the following way [J7]] : 

(a) = <pcos(2/x) , (tiq) = <psin(2/x) , where x is a spatial coordinate. (2.9) 

In the limit / — > the v.e.v. is a constant equal to <p, and the ordinary parity doublet model is 
realized. Inserting this Ansatz the Lagrangian for the baryons as well as for the mesons changes. 
For the potential Ym a further contribution 2f 2 (p 2 follows. The contribution arises from the kinetic 
term for the pions and sigma. As expected, this term suppresses the formation of a CDW. Due 
to the fact that the CDW is suppressed in the vacuum the parameters remain the same as in the 
case without CDW. Also the explicit symmetry breaking term get slightly modified at finite baryon 
chemical potential. 

For the baryons the modifications are more demanding. The Lagrangian has now a explicit 
coordinate space dependency. By transforming the fields \ffi and Xfc as [[P 



Yi -> Yi exp[-iftT 3 /x] , yi -»■ exp[-iy 5 z 3 fx)\i/i , (2-10) 
y/2 -»• ^ 2 exp[+;7 5 T 3 /x] , x\f 2 -> exp[+ij 5 r 3 fx] y 2 , (2-11) 

the Lagrangian takes the following form: 

££ B =\j/ii^¥i + ^17175^3/^1 + W$Y2 - ¥27i75^f¥2 

-^Si^ViYi-^h^ViYi + m^YsYi-YiYsVi) ■ (2-12) 

For illustrative purpose the additional COq dependency and baryon chemical potential jj, is not 
shown, but it is simple to generalize the expressions. The diagonalization of the Lagrangian leads to 
a splitting of the energy eigenstates. In the presence of the CDW four different energy eigenstates 
emerge, in contrary to the homogeneous case where only Af and Af* are present. The different states 
can be written down in the general form E^ip) = \j p 2 + mi c (pi) 2 , k = 1 . . .4, where m^pi) 2 is a 
momentum dependent mass which has to be calculated numerically. The grand canonical potential 
is given by: 

a/V = 2/V + l -X(p A - l -m 2 (p 2 -eq>- l -mlm 2 

4 2 r°° 1 

+ J JPi®[^* -E k (pi)]-[3^E k ( Pl ) 2 -(^i*) 3 -2E k ( Pl ) 3 ] , (2.13) 
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with jit* = jU — g a (0o- The behavior at nonzero baryon chemical potential can be found by extrem- 
izing the grand canonical potential with respect to the dynamical degrees of freedom: 

n ± <m/v) n ± d{si/v) nl d(a/v) 

dwo ' U " dcp ' df • 1 ' 

The parameters g m , A, m, e and m a are fixed using the known vacuum values for the masses 
of the pions and omega, as well the pion decay constant (2X(p 2 = m 2 — m\, 2m 2 = m 2 — 3m 2 , 
e = m 2 n f % ) and the conditions for nuclear matter saturation: 



1 

,1 ' ' fm 



— -m w = -16MeV, andp = 0.16 — . (2.15) 



The first equation ensures a binding energy per nucleon and the second the density at the minimum 
(it corresponds to the baryon chemical potential of jx = 923 MeV). The parameter mo remains as a 
input parameter. Using in addition the masses for N and N* the parameters gi and §2 can be fixed. 

3. Results 

The phase diagram of the parity double model shows a rich structure at dense hadronic matter. 



Since the model is based on Ref. [13] we use their parameters as a starting point. For different val- 
ues of the mass m^* and different values for the parameter mo the chiral condensate (p is calculated 
as a function of the baryon chemical potential pL. A feature of the model is to achieve nuclear sat- 
uration. For all combinations of the values for m^* and mo the parameters m and g a are tuned in a 
way to fullfill the conditions for nuclear saturation at a baryon chemical potential of pL = 923 MeV. 
The results are shown in Fig. 1. On the left hand side m^* is low, m^> = 1200 MeV and on the right 
hand side is high, m^* = 1500 MeV. For both cases the parameter mo is varied from 600 MeV to 
700 MeV and finally to 800 MeV. A general result for all parameter sets for the behavior of 9 is 
the presence of three different phases. 

For small /I the chiral symmetry is broken and the value for the condensate <p remains constant 
at the vacuum value f„. For an intermediate range for /I in the order of /I = 923 MeV one obtains 
a first order phase transition. The value for the condensate (p drops from the vacuum value f„ to 
30 — 60 MeV, where the ground state is still homogeneous and the CDW is not favored. In this 
range nuclear matter is formed. It shows that for the formation of homogeneous nuclear matter 
the mass of the pion plays a crucial role. Sending m n — > 0, the CDW seems to be always the 
favored ground state. This is unaffected by the choice of m^. For large /I a second first order phase 
transition occurs. The CDW is now always the favored ground state. Also for asymptotic large 
values for /i chiral symmetry will not be restored and the CDW remains the favored state of matter. 

For m^v* = 1500 MeV (left panel of Fig. 1) the second phase transition shows a strong depen- 
dence on the parameter mo. For small values of mo the transition happens at large pL and for a larger 
mo the intermediate homogeneous ground state shrinks to a small range of p.. The value for q> after 
the second phase transition also displays a dependence on mo, still the overall behavior is the same. 

For m;v* = 1200 MeV (right panel of Fig. 1) the situation is slightly different: the behavior of 
<p shows no mo dependency after the second phase transition. Furthermore for mo = 600 MeV and 
mo = 700 MeV the critical jj. is almost the same. 
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Figure 1: The chiral condensate <p as a function of the baryon chemical potential pi for different values of mo 
and niN* . On the left panel with a large mass for the chiral partner of the nucleon m^* = 1500 MeV and on 
the right panel a small mass mjv* = 1200 MeV. The red color indicates the ground state to be homogeneous, 
while the green that the CDW is favored. 
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Figure 2: The left plot shows the dispersion relation for the CDW. The spectrum shows four different energy 
levels in contrast to the homogeneous ones. The right shows the relative density of the states at finite ji. 



Another interesting point is the dispersion relation in the inhomogeneous phase, see Fig. 2 left 
panel. It shows the that the spectrum splits from two energy levels in the homogeneous phase to 
four in the inhomogeneous one. Moreover, for the lower two the energy decreases with increasing 
momentum and the upper two increase with increasing momentum. The parameters are the same 
as in the case for mo = 800 MeV and yi is slightly larger than the critical value for the transition to 
the CDW. Using the same parameters as in the left plot of Fig. 2 the relative densities are shown. 
Within the homogeneous region the transition to nuclear matter is found. 

4. Summary 

Inhomogeneous condensation is a relevant feature of the parity doublet model coupled to the 
linear sigma model with vector mesons. In a first step it has been shown that for m K = 139 MeV 
homogeneous nuclear matter is obtained. At a large baryon chemical potential the CDW is the 
favored state of matter. Sending the pion mass to zero (chiral limit) the picture changes. It is now 
no longer possible to find a homogeneous ground state. Furthermore, in both cases chiral symmetry 
does not get restored at high baryon chemical potential and the chiral condensate even increases 
for increasing density. 
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For the future the model has to be further extended. The CDW is just one of many different 
realizations of inhomogeneous condensation. A more general ansatz should be used [§]. The ex- 
tended linear sigma model (eLam) [|l^] has proven to be robust in describing vacuum phenomenol- 
ogy. Therefore combining both the parity doublet and the eLam would allow to constrain the ex- 
istents of CDW. Also the effect to the onset of the CDW regime of exotic matter like tetraquarks or 



glueballs [ 14 ] has to be tested in a more general framework. 
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